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ABSTRACT 

This paper discusses properties of the European Large Area ISO Survey 15 um quasars 
and tries to establish a robust method of quasar selection for future use within the 
Spitzer Wide-Area Infrared Extragalactic Survey (SWIRE) framework. The impor- 
tance of good quality ground-based optical data is stressed, both for the candidates 
selection and for the photometric redshifts estimates. Colour-colour plots and template 
fitting are used for these purposes. The properties of the 15 um quasars sample are 
studied, including variability and black hole masses and compared to the properties 
of other quasars that lie within the same fields but have no mid-infrared counterparts. 
The two subsamples do not present substantial differences and are believed to come 
from the same parent population. 
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1 INTRODUCTION 

The Europe an Large Area ISO Survey (ELAIS) 

iOliver et al.l [2000) is the largest survey performed 
with the Infrared Space Observatory (Ke ssler et alJ ITflflff) 
at 6.7, 15, 90 and 175 um and result ed in the delivery of the 
large st catalogue of any ISO surv ey <| Rowmi^^binson et alJ 
12004) from both the ISOCAM ICesarskv et al.lll996h and 
ISOPHOT llLemke et alJ Il996ft instruments. In partic- 
ular, the 15 um survey (performed with the ISOCAM 
instrument) covers an area of ~12 deg 2 , divided into four 
main fields (Nl, N2, N3, and SI) and several smaller 
areas. 15 u m observations in th e four main fields were 
analysed by I Vaccari et al.l ll2004t) . providing a catalogue 
of 1923 sources detected with S/ N > 5 over 10.85 deg 2 . 
The Final Band-Merged Catalogue llRowan- Robinson et alJ 
2004) combines all source lists at different wavelengths and 
redshifts obtained to date in ELAIS fields. This catalogue 
comprises a total of 3523 entries with about one third 
having spectroscopic identifications. 

Due to the fact that complete spectroscopic follow-up 
is usually not feasible over large and deep fields, one needs 
to use tools for detecting quasar candidates using photo- 



metric data only. As part of our study of mid-infrared (IR) 
quasars we present the results of two independent quasar 
candidates selection techniques, one based on colour-colour 
diagrams, and the other one template fitting, and try to 
improve the selection including IR constraints. In recent 
years, with the available multicolour surveys, quasar photo- 
metric redshift methods have been developed, yielding how- 
ever somewhat less reliable results than the ones for ga laxies 
(e.g. lHatziminaoglou et alJl20 00: Rich ards et alJ 12001'). For 
the purposes of this work, the template fitting technique 
(Hatz iminaoglou et a l. 2000) is applied on the two different 
data sets available for the studied fields (SDSS and WFS). 
All methods and results described in this paper can be di- 
rectly applied to the Spitzer Wide- Area I nfrared Extragalac- 
tic Survey ( SWIRE: ILonsdale et alJl2003ll . 

The layout of this paper is as follows. In Section [5] a 
brief description of the optical data available for this work 
is given. Section [3] deals with variability issues. In Section 
[I] a description of the selection of quasar candidates us- 
ing optical and IR properties is given, based on two differ- 
ent methods. Photometric redshifts for the spectroscopically 
confirmed quasars in the two data sets (SDSS and WFS) are 



2 Afonso-Luis et al. 



estimated and the results obtained for the two photometric 
systems are discussed. Section |^| compares the results for 
the sources with and without IR counterparts, in terms of 
their statistical properties and black hole (BH) masses. Our 
conclusions are presented in Section |S| 



2 THE 15 pm QUASAR SAMPLE AND 
RELATED OPTICAL DATA 

In the present work we study the type I quasars detected 
by ISO at 15 pm in two of the ELAIS fields, Nl and N2. 
Throughout this work, the traditional (but conservative) 
requirement of a quasar to be point-like has been set and 
sources classified as extended based on their r-band mor- 
phology are not taken into account. The morphological selec- 
tion is made in order to avoid contamination by galaxies of 
the quasar candidates samples discussed in Section0]and has 
no impact on the validity of the results presented here. Out 
of the 1056 sources in Nl and N2 contained in the ELAIS 
15 /m i Final Analysis Catalogue Version 1.0 llVaccari et alJ 
2004) , 849 sources wer e iden tified in INT WFS images by 
iGonzalez-Solares et al.l <|2004l). the non- identification being 
due either to incomplete optical coverage of the 15 pm fields 
or the optical limiting magnitude. 

The ELAIS Nl and N 2 have been fully covere d by the 
Wide Field Survey (WFS; iMcMahon et al] J200ll) . carried 
out with the prime focus Wide Field Camera (WFC) at 
the 2.5m Isaac Newton Telescope (INT) at La Palma. The 
survey consists of single 600s exposures in five bands U, 
g, r, i and z down to AB magnitude limits (5a limits for 
a point source) of 24.1, 24.8, 24.1, 23.6 and 22.4, respec- 
tively. The AB corrections for the conversion from the (orig- 
inal) Vega magnitud es to AB magnitudes have been com- 
puted using HyperZ feolzonella et all lioOO) and are: 0.751, 
-0.063, 0.165, 0.407 and 0.534, for U, g, r, i and z, respec- 
tively. Out of the 849 sources detected in 15 pm in Nl and 
N2, there are 110 point sources (SExtractor CLASS_STAR 
^ 0.9), excluding saturated s tars, with emission at 15 pm 
iGonzalez-Solares et alj|2004h . 

In addition to WFS, the Sloan Digital Sky Survey has 
valida ted and made publicly available its Data Release 1 
(DR1: lAbazaiian et al.l l 20 03). partially covering our fields. 
SDSS consists of five-band (u, g, r, i, z) imaging data cov- 
ering 2099 deg 2 , 186.240 spectra of galaxies, quasars, stars 
and calibrating blank sky patches selected over 1360 deg 2 
of this area, and catalogues of measured parameters from 
these data. The imaging data reach a limiting AB magni- 
tude dVanden Berk et al.ll2003h of r ~ 22.6 (95% complete- 
ness limit for stars). Among the 15 pm sources in the areas 
covered by the SDSS DR1 photometry, 82 have been mor- 
phologically classified as point sources (SDSS OBJC _TYPE 
= 6). In fact, the partial coverage of Nl and N2 provides 
photometry and spectra for a variety of sources. More par- 
ticularly, among the 36 SDSS spectroscopically confirmed 
quasars lying in the Nl and N2 fields, 16 have been de- 
tected at 15 pm. Note, however, that the areas covered by 
spectroscopy do not exactly coincide with those covered by 
photometry and are, in fact, smaller. Two morphologically 
extended low-redshift quasars (0.214 and 0.245) with r-band 
magnitudes of 18.25 and 18.58, have been excluded from the 
spectroscopically confirmed quasar sample, for the reasons 



1 

• SDSS 


OWYFFOS 


I I 

• 


- 


o 
• 


z < 
0.5 
1 .5 


0.5 

<7 < 1 5 


• 


- 


• • 


• 


* 

o 






oo 

s •. 


. o. 








' ° o 

• 


• 




I 





0.2 0.4 0.6 

r— band flux (mJy) 



Figure 2. r-band versus the 15 um flux for the quasar sam- 
ple. Filled (open) circles represent objects with SDSS (WYFFOS) 
spectra. Symbol size decreases with increasing redshift. 

mentioned at the beginning of this Section. Throughout this 
work and unless otherwise stated, all magnitudes will refer 
to the AB system. 

Fig.[I]shows the coverage of the ELAIS Nl and N2 fields 
by the WFS and the SDSS DR1 photometric and spectro- 
scopic data. 

Finally, as a part of the extensive program of ground 
based spectroscopy associated with the ELAIS survey, there 
have been follow-up observations of ELAIS 15 pm sources 
in the Nl and N2 fields (Perez-Fournon et al. in prepara- 
tion) using the WYFFOS multi-fibre spectrograph, on the 
William Herschel Telescope (WHT) at La Palma. This ob- 
serving run added nine previously unknown quasars to the 
spectroscopic sample in the regions of Nl and N2 not cov- 
ered by the spectroscopic DR1. 

Therefore, a total of 25 spectroscopically confirmed 
quasars have been identified among the final ELAIS 15 pm 
catalogue and their properties are described in Table Q Fig. 
[5] shows the 15 pm flux versus r-band flux for the quasar 
sample. 



3 VARIABILITY 

Variations in the luminosity of quasars have been observed 
from X-ray to radio wavelengths, with timescales of min- 
utes to years. The majority of QSOs have continuum vari- 
ability on the order of 10% on timescales of months to 
J | II i 

years (Vandcn Berk et al. 2003). Furthermore, recent obser- 
vations of radio-quiet quasars indicate that more than 80% 
show long-term (month t o year) variability w ith amplitudes 
up to half a magnitude l|Huber et al.l 120021 . Variability is 
wavelength dependent. The continuum {F v oc v~ a ) tends to 
get harder (the spectral index, a, decreases) as the quasar 
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Figure 1. Coverage of the ELAIS Nl and N2 fields (small inner rectangles) by the WFS and the SDSS DR1 photometric and spectroscopic 
data. 



Table 1. ELAIS 15 um spectroscopically confirmed quasars, r-band magnitude in column 5 is the SDSS r-band magnitude apart from 
the objects marked with an asterisk. These objects are outside the spectroscopic coverage of D R1 and th e magnitude shown is the 
WFS r magnitude. The last column (reference) refers to the origin of the spectra ( 1: lAbazaiian et alJl20oH: 2: Perez-Fournon et al. in 
preparation) 



ELAIS ID 


RA (optical) 


Dec (optical) 




z 


r mag 


15 pm flux (mjy) 


z Reference 


ELAISC15 


.J160250.9+545057 


240.71234131 


54.84947205 


1 


1971 


19.22 


2.1820 


SDSS 1 


ELAISC15 


J160522.9+545613 


241.34640503 


54.93708420 





5722 


18.94 


2.2160 


SDSS 1 


ELAISC15 


J160623.5+540556 


241.59822083 


54.09888840 





8766 


17.62 


4.2930 


SDSS 1 


ELAISC15 


J160630.5+542007 


241.62754822 


54.33544159 





8205 


18.73 


3.4460 


SDSS 1 


ELAISC15 


J160638. 0+535009 


241.65779114 


53.83572006 


2 


9426 


19.77 


1.6010 


SDSS 1 


ELAISC15 


J161007.2+535814 


242.52960205 


53.97055435 


2 


0317 


18.86 


1.7400 


SDSS 1 


ELAISC15 


J163702. 2+413022 


249.25930786 


41.50616837 


1 


1783 


19.11 


2.0920 


SDSS 1 


ELAISC15 


.J163709.2+414031 


249.28890991 


41.67519760 





7602 


17.20 


8.4090 


SDSS 1 


ELAISC15 


J163739.3+414348 


249.41436768 


41.72999954 


1 


4136 


18.94 


1.0610 


SDSS 1 


ELAISC15 


.J163847.5+421141 


249.69760132 


42.19494629 


1 


7786 


18.93 


1.7420 


SDSS 1 


ELAISC15 


J163915.9+412834 


249.81591797 


41.47602844 





6919 


19.05 


3.3730 


SDSS 1 


ELAISC15 


J163930.8+410013 


249.87844849 


41.00380707 


1 


0515 


18.23 


1.2580 


SDSS 1 


ELAISC15 


.J163952.9+410346 


249.97023010 


41.06244278 


1 


6050 


18.58 


2.0940 


SDSS 1 


ELAISC15 


J164010.1+410521 


250.04225159 


41.08955383 


1 


0990 


17.01 


9.5570 


SDSS 1 


ELAISC15 


J164016.0+412102 


250.06701660 


41.35038757 


1 


7570 


18.44 


2.0900 


SDSS 1 


ELAISC15 


J164018.4+405812 


250.07641602 


40.97030640 


1 


3175 


18.13 


3.7700 


SDSS 1 


ELAISC15 


.J161521.8+543148 


243.84077454 


54.53016663 





4737 


18.24* 


3.492 


WYFFOS 2 


ELAISC15 


.J161526. 7+543004 


243.86094666 


54.50175095 


1 


3670 


19.35* 


1.457 


WYFFOS 2 


ELAISC15 


.J161543. 5+544828 


243.93133545 


54.80799866 


1 


6920 


18.22* 


2.107 


WYFFOS 2 


ELAISC15 


.J163425. 2+404152 


248.60472107 


40.69794464 


1 


6840 


18.37 


3.1690 


WYFFOS 2 


ELAISC15 


J163502. 7+412953 


248.76176453 


41.49808502 





4727 


18.08 


2.1730 


WYFFOS 2 


ELAISC15 


J163531. 1+410025 


248.87956238 


41.00761032 


1 


1500 


18.81 


1.3990 


WYFFOS 2 


ELAISC15 


.J163533. 9+404025 


248.89176941 


40.67380524 





5340 


19.72 


2.7300 


WYFFOS 2 


ELAISC15 


J163553.5+412054 


248.97354126 


41.34883118 


1 


1950 


19.39 


2.7260 


WYFFOS 2 


ELAISC15 


J163634.4+412742 


249.14343262 


41.46200180 





1711 


18.50 


4.1540 


WYFFOS 2 



gets brighter, which means that the variations are larger at 
shorter wavelengths. 

The SDSS imaging strategy consists in observing in al- 
most simultaneous mode the same field in the five differ- 
ent bands. The imaging strategy of the WFS was, however, 
different, and the same fields were observed in the differ- 
ent filters in timescales ranging from a few months to more 
than a year. Fig. |2l shows the differences of the magnitudes 
for the point sources lying in the common area in the two 



photometric systems. The magnitude differences show very 
small dispersions down to a certain magnitude for almost all 
point sources that are identified as stars but large variations 
appear in the case of many of the confirmed quasars (red 
crosses) in all filters but especially in the U- and g-bands, 
indicating variability related issues. 

It might be argued that the variations are due to the 
different photometric systems, however close they may be. 
For this purpose, a simulated quasar catalogue was created, 
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Figure 3. Magnitudes differences between SDSS and WFS, the black points represent the point sources, mostly stars, while the red 
crosses indicate the 73 confirmed quasars. The photometric errors for the quasars are lower than 0.05 magnitudes in all bands and the 
dispersion in their magnitudes is evident. 



based on the colours of the SD SS composite quasar spec- 
trum JVanden Berk et al-lEoOll) . including all 10 filters and 
for redshifts spanning from to 6 and the magnitudes were 
compared. Fig. |3] illustrates the expected magnitude differ- 
ences as a function of magnitude for redshifts in the interval 
[0,4], where all our spectroscopically confirmed quasars lie, 
and for the filters U, g and z. Filters r and i were omit- 
ted for clarity as the points largely overlap with the rest 
of the points. Even though magnitude variations are pre- 
dicted, their amplitude is significantly smaller than the ones 
observed. We therefore conclude that not all observed vari- 
ations between the SDSS and WFS photometry can be ex- 
plained by the differences of the photometric systems and 
that some must be a result of quasar variability. 

This strongly suggests that quasar candidates selection 
and photometric redshift estimates are likely not to be ac- 
curate for all objects when WFS or other similar data are 
used and additional methods have to be looked for. 



the quasar candidates. The point source template library 
comprises quasar and stellar templates and the observed 
SED of each object is compared to the one computed by 
convolving each t emplate with the filters tra nsmissions (for 
further details see lHatziminaoglou et all200d and 2002). Ml 
is expected to give higher confirmation rates (i.e. number of 
real quasars over then number of quasar candidates) at low 
and intermediate redshifts but its efficiency greatly depends 
on the photometric system used and can not be applied as 
such when other filters are used. M2 is subject to higher 
contamination from sources other than quasars (e.g. stars) 
but is expected to have a much better efficiency at high red- 
shift and can be used for any filter combination. In order to 
improve the results, we also make use of the IR information 
available for the 15 um sources keeping in mind that the 
combination of all these techniques will be applied in the 
near future in the framework of SWIRE. 



4 QUASAR CANDIDATES SELECTION 

For the purpose of identifying high probability quasar can- 
didates from the available datasets we used two techniques, 
a combination of colour-colour diagrams (hereafter method 
1 - Ml) and the t e mplat e fitting method (MS). Ml, based 
on lRichards et aT] J2Q02), consists of a colour-colour selec- 
tion algorithm trained using the SDSS Early Data Release 
JStouehton et al.ll2002t) . for low- to intermediate redshift 
(~ 2.5) objects. M2 is the standard template fitting that 
simultaneously provides a photometric redshift estimate for 



A first test is made using the more reliable SDSS pho- 
tometry fSection l4.1|l but an attempt of selecting candidates 
based on the WFS photometry will also be presented (Sec- 
tion ^^J- O ur t es t sample consists of 21 spectroscopically 
confirmed quasars with available SDSS photometry and 25 
with WFS photometry, which include all 21 from SDSS (see 
Table Fine-tuning the method for WFS data is very im- 
portant as large part of the SWIRE fields have been observed 
by the WFS. 
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Figure 4. Theoretical magnitude differences for the filters U, g, 
and z of the SDSS and WFS photometric systems, in the redshift 
interval [0,4], Filters r and i have been omitted for clarity. 



Table 2. Comparison of the completeness and confirmation rate 
yielded by Ml and M2 on a sample of 82 SDSS point sources with 
f 5 pm emission. The second column shows the number of quasar 
candidates (N c ) and the number of spectroscopically confirmed 
quasars among them (Nf). 



Method 


N c (N f ) 


Completeness 


Conf. Rate 


Ml 


26 (19) 


90% 


73% 


M2 


33 (19) 


90% 


58% 


M2&. CI 


30 (19) 


90% 


63% 



4.1 Quasar Candidates with SDSS Photometry 

A sample of 82 15 pm ELAIS sources identified as point 
sources from their optical photometry (SDSS OBJC.TYPE 
= 6) with r-band magnitudes brighter that 22.6 has been 
selected from the SDSS photometric catalogue. The optical 
magnitude cut has been imposed in order to avoid spurious 
detections and large photometric errors. 

The 15 pm information can be used in order to im- 
pose IR conditions in the selection of quasar candidates. 
Stars, galaxies, and AGN all have different optical to mid-IR 
slopes, with stars typically having lar ger optical than mid-IR 
fluxes (jGonzalez-Solares et al.ll2004). Furthermore, accord - 
ing to models of galaxies in the IR jRowan-Robinsonll200ll) , 
quasar mid-IR fluxes are some 10 to 100 times larger than 
their optical ones. Taking this into account, one can impose 
additional constraints on the selection criteria requiring a 
mid-IR to optical flux ratio of at least 10 for quasar can- 
didates (hereafter condition CI). This condition allows the 
removal of three quasar candidates selected by M2 that have 
mid-IR to optical fluxes those of stars. Fig. illustrates the 



positions of the different objects types and the regions where 
quasars and quasar candidates selected by the two methods 
lie. 

In order to distinguish between quasars and galax- 
ies, one can make a combined use of optical and opti- 
cal/IR colours taking advantage of the fact that quasars 
up to a redshift of ~ 3 are typically bluer than galaxies 
jGonzalez-Solares et al]l2004l) . Furthermore, 97.55% of all 
spectro scopically confirmed q uasars in the DR1 quasar cat- 
alogue dSchneider et alj |2003'l have r — i < 0.52 (hereafter 
C2). Fig. [(J shows the distribution of stars (lower mid-IR to 
optical fluxes), galaxies and quasars (bluer than galaxies in 
general). Marked in red are the spectroscopically confirmed 
quasars (left panel) and the candidates selected by meth- 
ods Ml and M2 (middle and right panels, respectively). As 
can be seen, all spectroscopically confirmed quasars form a 
clump bluewards of r — i ~ 0.52. In this particular case, C2 
does not improve the results of any of the methods but will 
be used further on. 

If N c is the number of quasar candidates stemming 
from the identification technique, Nf the number of real 
quasars among the candidates, and N e the number of 
expected (based on models) or known (based on com- 
plete observations) quasars, completeness and confirma- 
tion r ate can be defined as NjjN e and Nf/N c , respec- 
tively jHatziminaoglou et al] 120001 . Table H compares the 
two quantities yielded by the two methods. Both methods 
give the same completeness but the colour-colour selection 
favours the confirmation rate. Note, however, that the val- 
ues given here for confirmation rate are lower limits. A sub- 
stantial number of candidates (some 30%) are fainter than 
the SDSS spectroscopic completeness limit or lie outside the 
area covered by spectroscopy and therefore their nature is 
unknown. For the same reason, the value of the complete- 
ness is also indicative as a complete spectroscopic coverage 
would alter both Nf and N e . 



4.2 Quasar Candidates with WFS Photometry 

From the objects morphologically identified as point sources 
(SExtractor CLASS_STAR ^ 0.9) with detections at 15 pm 
110 obj ects are part of the ELAIS Fi nal Band- merged Cat- 
alogue iRowan- Robinson et aljl2004t) . Ml yields 27 quasar 
candidates with 18 spectroscopically confirmed, while M2 
finds 63 candidates with 25 spectroscopically confirmed (Ta- 
ble 0. The two candidate lists have 27 objects in common, 
18 of which are spectroscopically confirmed quasars. Tak- 
ing into account condition CI, 7 out of the 63 candidates 
proposed by M2 can be safely discarded as stars. Applying 
C2 on the remaining 56 candidates, we discard another five 
as more likely to be galaxies. Fig. |7] shows the positions of 
confirmed quasars and candidates, similar to Fig.|Bj but for 
the WFS photometry. 

The results obtained using WFS photometry are sum- 
marised in Table [3] As already mentioned, the values of 
completeness are underestimated due to lack of complete 
spectroscopic coverage. 

After thorough consideration we reach the conclusion 
that WFS data can be used in order to reliably obtain quasar 
candidates despite the issues raised by variability, especially 
if one combines Ml and M2 with the constraints imposed 
by the objects' IR properties. 
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Figure 5. 15 um flux (mjy) vs SDSS r-band magnitude for the point sources (stars) within the parts of ELAIS Nl and N2 areas covered 
by SDSS. The red squares indicate the confirmed quasars (left panel), Ml quasar candidates (central panel) and M2 quasar candidates 
(right panel). The blue lines represents constant 15 um to r-band flux, with values 0.1, 1., and 10. 




Figure 6. 15 um to r-band flux vs SDSS r — i for the point sources (stars) and extended sources (open circles). Red squares indicate 
spectroscopically confirmed quasars (left panel), Ml quasar candidates (central panel) and M2 quasar candidates (right panel). 
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Figure 7. 15 um to r-band flux vs WFS r — i for the point sources (stars) and extended sources (open circles). In red squares, the 
spectroscopically confirmed quasars (left), Ml quasar candidates (centre) and M2 quasar candidates (right). 
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Table 3. Comparison of the completeness and confirmation rate 
yielded by Ml and M2 on a sample of 110 WFS point sources with 
15 um emission. The second column shows the number of quasar 
candidates (N c ) and the number of spectroscopically confirmed 
quasars among them (Nr). 



Method 


N c (JV» 


Completeness 


Conf. Rate 


Ml 


27 (18) 


72% 


67% 


M2 


63 (25) 


100% 


42% 


M2U CI 


56 (25) 


100% 


45% 


M2h CI & 


5 C2 51 (25) 


100% 


49% 



4.3 Quasar Photometric Redshifts 

For estimating the photometric redshifts we applied a stan- 
dard template fitting procedure using synthetic quasar spec- 
tra consisting of a power law co ntinuum and emission lines o f 
fixed equivalent width values dHatziminaoglou et alJ boOO). 
on the sample of 73 spectroscopically confirmed quasars with 
both SDSS and WFS photometry available. The results are 
rather reliable when using SDSS photometry (53 out of 73 
Az ^ 0.2, i.e. 73% ), but get worse when using the WFS 
photometry (29 out of 73 for Az < 0.2, i.e. 30% ). Note, how- 
ever, that all seven objects with spectroscopic redshifts lower 
than 0.3 have been assigned the wrong photometric red- 
shifts. These objects, however, are extended (OBJCLTYPE 
= 3) and, therefore, their magnitudes (psf magnitudes for 
SDSS and core magnitudes for WFS) must be contaminated 
by the light of the host galaxy. If we consider only the objects 
with spectroscopic redshifts higher than 0.3, the numbers of 
good identifications become 80% and 44%, for SDSS and 
WFS, respectively. Furthermore, all 30 objects that were 
assigned correct photometric redshifts using the WFS pho- 
tometry have also correct photometric redshifts when SDSS 
photometry is used. Fig. |H| illustrates the results obtained 
for SDSS (left panel) and WFS (right panel) photometry. 

The fact that the bands for WFS photometry were 
taken in different periods of time (especially the [/-band, 
which was taken more than two years apart in some cases) 
lead us to the conclusion that variability might be the basis 
of the discrepancy problem we encountered, as described in 
Section [3| 

A case-by-case study of all objects that have correct 
photometric redshift estimates with SDSS photometry and 
bad estimates (Az > 0.25) with WFS photometry showed 
that they all lie in the redshift range [0.3,2.0] and their dif- 
ferences in magnitudes span a much larger range than the 
one expected due to the filters' differences. One finally con- 
cludes that the wrong photometric redshift assignments are 
most probably due to variability, as argued earlier. 



5 COMPARISON BETWEEN IR DETECTED 
AND NON-DETECTED QUASARS 

36 spectroscopically confirmed quasars lie in the parts of the 
ELAIS Nl and N2 fields covered by the SDSS spectroscopic 
data release, with 16 of them detected at 15 um, as already 
mentioned. Considering the h omogeneous way th e SDSS 
quasar candidates are selected dRichards et al . 2002) we can 



ask is if the two subsamples (16 IR and 20 non-IR emitters) 
have the same properties or if the IR detected quasars are, in 
some way, different. The colour-redshift (Fig|5J and colour- 
colour (Fig llOl diagrams of the two subsamples do not show 
any differences. 

A two-sided Kolmogorov-Smirnov test on the BH mass 
distributions of the two subsamples (seen in the right panel 
of Fig, 1111 resulted in a larger than 90% probability for them 
to come from the same population. The same test gave a 35% 
probability for the redshift distributions to be representa- 
tive of from the same population. The largest deviation is 
noticed in the r-band magnitude distributions (left panel of 
Fig. 1111 . The objects not detected by ISO at 15 um are on 
average half a magnitude fainter in the r-band. Fig. ^indi- 
cates a possible correlation between the 15 pm fluxes of the 
SDSS quasars and their optical (r-band) fluxes, suggesting 
that the lack of 15 um counterparts could be due to their 
fainter magnitudes. Deeper IR observation would, presum- 
ably, provide counterparts for the remaining objects, since 
for a given covering fraction and varying bolometric lumi- 
nosity, brighter optical AGN should have a higher mid-IR 
emission. 



5.1 Quasar Black Hole Masses 

The principal assumption underlying the BH virial mass es- 
timate is that the dynamics of the Broad Line Region (BLR) 
are dominated by the gravity of the central supermassive 
black hole. Under this assumption an estimate of the cen- 
tral BH mass, Mbh, can be Mbh — RblrV 2 /G; where 
Rblr is the radius of the BLR and V is the velocity of the 
line-emitting gas, traditional ly estimated from the FWHM 
of the Hp emission line (see lKaspi et al.ll2000l) . For quasars 
with redshifts hi gher than typically z ~ 0.8, when Hf3 is no 
longer available. iMcLure fc JarvisTteOOSl) suggested the use 
of Mgll as an estimator of the B H mass. More precisely, th e 
BH mass is computed as follows jMcLure fc Dunlodl2004h : 

Mbh _ ? / AL 5100 \ - 61 / FWHM{H P ) V 

M Q ■ U0 37 W7 ^ km/ sec J ( ' 

and 

M bh _ / AL 3 ooo \ - 62 / FWHM{MgII) \ 2 

M @ Vl0 37 WV ^ km/ sec ) 1 ' 

For the sample of 36 SDSS quasars in Nl and N2, there- 
fore, low redshift (up to ~ 0.8) mass estimates are based on 
the correlation between the H/3 and the monochromatic lu- 
minosity at 5100 A (Eq. while for higher redshift ones 
the correlation between Mgll and 3000 A luminosity (Eq. 
|2J is used. For objects with redshifts higher than 2, all es- 
timators fall outside the covered wavelength range and BH 
masses can no longer be computed. Only objects with SDSS 
spectra were used for the BH mass computation as line and 
continuum of all SDSS spectra are measured in a consistent 
way by the SDSS pipeline. 

Fig. ll2l shows the distribution of BH mass as a function 
of redshift. Red marks the objects for which H/3 was used 
(with redshifts typically lower than ~ 0.8), while black illus- 
trates the ones for which Mgll was considered. Filled circles 
denote the objects with 15 pm emission while triangles corre- 
spond to SDSS objects within the ELAIS fields but without 
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Figure 8. Photometric versus spectroscopic rcdshifts for the 73 SDSS quasars within the ELAIS Nl and N2 areas, using SDSS (left) 
and WFS (right) photometry. The larger scatter in the right plot is most likely due to quasar variability, see text for details. 









1 




1.5 
t.O 
0.5 
0.0 



-0.5 





2 3 4 
Redshift 



1 .2 3 4 
Redshift 



12 3 4 
Redshift 



1 



Figure 9. Colour - redshift diagrams for SDSS DR1 quasars (dots), non-ELAIS quasars (triangles) and ELAIS quasar (filled circles). 
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Figure 10. Colour - colour diagrams for SDSS DM quasars (dots), non-ELAIS quasars (triangles) and ELAIS quasar (filled circles). 
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Figure 11. Histograms of r- band magnitude (left panel), redshift (middle panel) and black hole mass (right panel) for the ELAIS (black 
solid line) and non-ELAIS (red dotted line) quasars. 
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Figure 12. Black hole mass as a function of redshift for the 
ELAIS quasars with SDSS spectra (filled circles) and for the SDSS 
quasars without ELAIS counterparts (triangles). 



IR detections. BH masses show no differences between the 
two subsamples and follow the distribu tions found for the 
entire SDSS DR1 quasar catalogue by iMcLure fc Duniopl 
( 2004). The right panel of Fig. ITT1 shows the BH mass his- 
togram for the two subsamples, clearly indicating that both 
of them belong to the same parent population. 

For the objects for which the 15 um fluxes were avail- 
able, their IR luminosity was computed assuming an (fi, A) 
= (0.3, 0.7) cosmology. Fig. ll3l shows the 15 pm IR luminos- 
ity of these objects as a function of redshift (left panel) and 
as a function of the BH mass (right panel). Clearly, high- 
mass BH tend to produce higher IR luminosities. Further- 
more, the IR luminosity versus redshift distribution com- 



pares nicely to the findings of lPozzi et al.l i2003T) , indicating 
that type I AGN in the redshift range [ 1,21 have IR lum i- 
nosities higher than 10 12 L Q (see Fig. 6 of lPozzi et ai]l2003l) . 
proving thus that ISO only skimmed the brighter end of the 
luminosity function. 



6 DISCUSSION 

This paper discusses some properties of the ELAIS 15 pm 
quasars and tries to establish a robust method of quasar se- 
lection for future use within the SWIRE framework. The im- 
portance of good quality ground-based data is stressed, both 
for the candidates selection and for the photometric redshifts 
estimates. Colour-colour plots and template fitting are used 
for these purposes. Colour-colour plots give a higher confir- 
mation rate, as expected since they were trained on the same 
photometric system, while template fitting has the advan- 
tage of being independent of the filters used. The SDSS data 
have proved to be a reliable dataset, while the WFS data are 
somewhat less efficient due to the large timescales over which 
the imaging has been carried out and the variability issues 
this introduced. In order to correct for these effects, the IR 
and optical to IR properties of the objects have been taken 
into account imposing additional constraints on the quasar 
candidate selection techniques. The two subsamples of spec- 
troscopically confirmed quasars detected and non-detected 
at 15 um have been examined and they have shown no intrin- 
sic physical differences. Their non-detection is most proba- 
bly due to their fainter magnitudes, probably correlated to 
their higher redshifts. 

This work can be seen as a validation of the tools and 
methods that will be used in the framework of SWIRE or 
other similar IR surveys supported by ground-based optical 
data. SWIRE will survey six high-latitude fields (including 
ELAIS Nl and N2), totaling 50 square degrees in all seven 
Spitzer bands. One of the key scientific goals of SWIRE is 
to determin e the evolution of qu asars in the redshift range 
0.5^ z llLonsdale et al.ll2004l) . As we suggested in Sec- 
tion |5J with deeper IR observations one could detect all the 
quasars in the discussed fields down to the optical magnitude 
limit. In particular, the SWIRE 5a photometric sensitivity 
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Figure 13. 15 pm IR luminosity of these objects as a function of redshift (left panel) and as a function of the BH mass (right panel). 



is 0.0037 mjy in the IRAC 8 pm band and 0.15 mjy in 
the MIPS 24 pm band, much deeper t han the characteristi c 
depth of ~ 1 mjy of ELAIS 15 pm l|Vaccari et alJl2004) . 
Last but not least, the multitude of the SWIRE bands will 
allow for a much better galaxy-AGN separation in the mid- 
IR colour space, possibly making the optical morphological 

preselection of qu asar candidates obsolete. 

Models (e.g. lOranato fc j3anggej_ ll994l: iNenkova et all 
12002?) and observations (e.g. lElvis et alJll994l) suggest that 
quasar IR spectra are more or less flat (in vL v ) from ~ 1 
pm down to at least 25 pm, even in the absence of starburst 
activity. Therefore, the deeper and better quality SWIRE 8 
pm and 24 pm observations will allow the detection of high 
numbers of quasars and the easy adaptation of the tools 
presented here. The stellar contamination will probably be 
higher at 8 pm than at 15 pm but this problem will most 
likely be solved due to the multitude of IR bands, that will 
allow an easier separation of the galactic and extragalactic 
populations. 

As a last remark, we would like to stress that a robust 
candidate selection technique and subsequent photometric 
redshift estimates such as the ones presented here will be in- 
creasingly required by all future large area surveys, as spec- 
troscopic coverage will never reach the same completeness, 
neither in area nor in depth. 
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